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1. Introduction 

The dynamical evolution of ultracold atomic quantum gases driven far out of equilibrium is 
a subject of intensively growing interest. Precision measurement techniques for many-body 
observables have been and are being developed with vigourous effort. This technology has 
triggered a demand for progress in theoretically describing non-equilibrium quantum many- 
body dynamics of strongly interacting systems beyond mean-field approximations. Recent 
highlights of this development include, e.g., the variation and enhancement of the atomic 
interactions on the basis of Feshbach scattering resonances 0], the achievement of strongly 
correlated regimes within optical lattices |l2l, as well as the quench dynamics of spinor Bose- 
Einstein condensates Q. 

To obtain the full dynamics of a quantum many-body system, it is convenient to formulate 
an initial-value problem that describes the time evolution by means of in general coupled 
equations of motion for time dependent correlation functions and specific values for these 
functions at initial time. Alternatively, the time-evolution of many-body systems is very often 
described in terms of kinetic or transport equations JH. In general, the aim of quantum kinetic 
theory is to find evolution equations for distribution functions /(x, t) or /(p, t), interpreted, 
e.g., as particle number densities in x- or p-space. These distribution functions can then be 
used to derive various transport properties of the many-particle system as, e.g., current of 
charge or energy. For that reason, the evolution equations of / are generally referred to as 
transport equations. 

Kinetic descriptions usually neglect the effect of correlations between different times of 
the evolution, i.e., they build, to a certain extent, on a Markovian approximation. In particular, 
they neglect the initial dynamics directly after a change in the boundary conditions which 
drive the system out of equilibrium. This shortcoming is cured in dynamical approaches in 
which coupled equations of motion for the correlation functions are derived to describe the 
time evolution starting from a specific initial state. The built-up of correlations beyond the 
kinetic approximation, in these equations, is taken into account by means of non-Markovian 
integrations over the evolution history of the correlation functions. 

The focus of this work is set on the question to what extent transport theory can be 
used to quantitatively describe the dynamics of an ultracold Bose gas at different stages of 
its evolution from the initial to the equilibrium state. Determining the range of validity of 
transport equations by comparing to the full quantum dynamic evolution helps forming a solid 
foundation for the application of transport equations which, for practical problems, often is 
technically less demanding. 

In the present work we use the two-particle irreducible (2PI) effective action ^ which 
is derived from a functional integral representation of the quantum many-body system. From 
this action, self-consistent quantum-dynamic equations for the mean field as well as the two- 
point correlation functions are obtained |l6ll3. These have been applied successfully for the 
description of far-from-equilibrium dynamics and thermalization [[8l[6ll9l [T0l[TTl[T2l[T3l[T4ll . 
For our investigations, we numerically compare the time evolution determined with dynamic 
equations derived from the 2PI effective action with the corresponding transport equations. 
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These transport equations represent an approximation to the full dynamic equations. The 
derivation of transport equations has been discussed extensively in the literature. In the context 
of non-relativistic systems, in particular cold atomic gases cf., e.g., Refs. Il4l [TSl [T6l [TTl [TSl 
[H |20l im m [231 m. see also Refs. [[251 [23 |27l |28l IIS [30l [3ll El in the context of 
relativistic physics. A comparison of dynamic equations with their kinetic approximation 
for relativistic dynamics has been given in Refs. [[25l [33l [34l [TTl [3T1 [32l . Here, we derive 
transport equations in leading order (LO) and next-to-leading order (NLO) of a gradient 
expansion. This is essentially achieved in four steps: (a) the previously obtained 2PI dynamic 
equations for the two-point correlation functions which are expressed with respect to absolute 
coordinates are rewritten with respect to center and relative coordinates. This allows (b) for a 
gradient expansion with respect to center coordinates, (c) Furthermore, knowledge regarding 
the details of the initial state is neglected, (d) These steps allow for a Wigner transformation, 
i.e., a Fourier transformation with respect to relative coordinates, which yields the desired 
transport equations. We then use the results obtained from solving the 2PI dynamic equations 
to study the range of validity of the transport equations. 

While the solution of an initial- value problem requires to formulate a dynamical instead 
of a kinetic theory, a further important aspect is to properly take into account higher 
order correlations in the description of strongly correlated systems. Far from equilibrium, 
perturbative mean-field approximations such as in Gross-Pitaevskii |[35l and Hartree-Fock- 
Bogoliubov |[36ll theory cease to be valid descriptions. The dynamics of a self interacting 
ultracold gas is well-described by such mean-field approximations only as long as fluctuations 
and the interaction strength are small such that expansions in a small coupling parameter are 
justified. This is generally fulfilled for weakly interacting systems close to equilibrium. For 
strongly fluctuating systems this approach, however, becomes insufficient when requiring 
summations of infinite series of perturbative processes. But even close to equilibrium, 
fluctuations can play an important role on large time scales. 

In this article we consider a non-perturbative approximation which reaches substantially 
beyond the Hartree-Fock-Bogoliubov mean field theory. This non-perturbative approach is 
based on a systematic expansion of the 2PI effective action in powers of the inverse number 
of field components A/^ [[H [3. It has recently been used to describe the dynamics of an 
ultracold atomic Bose gas far from equilibrium |[T3l [T4| and its non-perturbative character 
has been reconfirmed in the framework of a functional renormalization group approach |[37l . 
The 2PI l/M expansion to next-to-leading order yields dynamic equations which contain 
direct scattering, memory and "off-shell" effects. It allows to describe far-from-equilibrium 
dynamics as well as the late-time approach to quantum thermal equilibrium. Recently, these 
methods have allowed important progress in describing the dynamics of strongly interacting 
relativistic systems far from thermal equilibrium for bosonic |[6l [38l [39l [H [40l [4T]| as well as 
fermionic degrees of freedom [[T0l[42l . 

Our article is organized as follows: In Section [2] we recall the functional description of 
the quantum many-body dynamics on the basis of the 2PI effective action and present the 
set of coupled dynamic equations for the two-point correlation functions. In Section [3] we 
derive, by transforming to Wigner coordinates and performing a gradient expansion, the set 
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of transport equations which we compare, in Section HI numerically with the 2PI dynamic 
equations. Our conclusions are drawn in Section [5l 



2. Dynamical evolution equations 

The evolution of single-particle observables such as the spatial density and momentum 
distributions are most conveniently described through the spectral and statistical two-point 
functions, denoted as p and F, respectively ll43l : 



Fab{x,y) = i({$„(x),<l>,(i/)})e= i({$a(x),$fe(y)}) (1) 

Pa,{x, y) = i([$a(x), = i([$a(x), (2) 

Here, (...) implies a trace over the density matrix describing the full quantum state at the 
initial time, and (. . .)c is a shorthand for the corresponding connected correlation function 
or cumulant. As defined in more detail below, we will consider a non-relativistic theory for 
a complex scalar field with J\f = 2 components describing an ultracold Bose gas, and we 
assume pointlike interactions between the atoms. Hence, all fields and n-point functions 
acquire indices a, b, ... G {1, 2}. 

From the above definition it is clear that the full propagator function G, defined as the 
expectation value of the time-ordered product of field operators, can be written in terms of the 
statistical and spectral functions as follows: 

Gab{x,y) = {Tc<l>a{x),<l>b{y))c 

= F{x,y) - ^p{x,y) sgn^ {xq - yo) . (3) 

Our aim is to solve an initial-value problem for the two-point function, such that we need to 
work with the Sch winger- Keldysh closed time path (CTP) C II44I1 from the initial time along 
the positive time axis to the largest time appearing in the product of field operators, and back 
to the initial time. Tc denotes time ordering along C, and the sign function sgn^^ evaluates to 1 
(— 1) for xo later (earlier) than yo along the CTP. 

While the spectral function p encodes the excitation spectrum of the theory, the statistical 
propagator F contains information about the occupation numbers. Roughly speaking, p 
makes explicit what states are available and F how often they are occupied. Far away from 
equilibrium, F and p are in general independent functions, while close to thermal equilibrium, 
they are connected by a fluctuation dissipation relation [|4ll45l. 

For Xq = yo the equal time commutation relations for the field operators $ imply 

= -ic^25(x-y). (4) 

The exact dynamical evolution equations for the statistical and spectral functions can be 
derived from the 2PI effective action [|6l |71 and read, for a vanishing field expectation value 

<Pa{x) = (Mx)) = 0: 



d,,Fix,y)= i(T2|[i^iB(x) + sW 



[X 



F{x,y) 
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+ / T.P{x,z)F{z,y)- T.^{x,z)p{z,y)\, (5) 

J zfi J zfl ) 

d,,p{x,y) = ia2|[iJiB(x) + sW(x)lp(x,y)+ [ ^^^x, z)p{z,y)] . (6) 



z,yo 



Here, a2 is the Pauli 2-matrix, and we have suppressed the explicit notion of field indices, i.e., 
all functions are understood to be 2 x 2-matrices. We use the notation = dzo [ d^z. 
Equations ^ and ^ represent an initial value problem, i.e., they give the time evolution of 
the functions F and p for given initial values F(to, x; to, y) and p(to, x; to? y)- The memory 
integrals on the right hand sides take into account the build up of higher-order correlations 
during the evolution and ensure the equations to respect causality. Since the equations are 
exact, they are equivalent to any other representation such as dynamical Kadanoff-Baym or 
Schwinger-Dyson equations for the time-dependent two-point functions. 

The spectral (T.^) and statistical (S^) parts of the self-energy, as well as the local energy 
shift are obtained from the proper self-energy S which sums all one-particle irreducible 
(IPI) diagrams, in a similar way as in Eq. ([3]): 

S(x,i/)= -iSW(x)5(x-y) + S(x,y), (7) 

y) = y) - ^^^{x, y) sgnc(xo - yo). (8) 

To close the set of dynamic equations © and ® , the self energy has to be specified. We 
derive the self energy from the beyond-one-loop contribution r2 to the 2PI effective action F, 

r[0, G] = r°— i°°p[0, G] + r2[0, g], (9) 

pone-ioopj^^ G] = S[<i)\ + ^Tr{ln(G-i)} + ]p:i{G^^[<i)]G} + const., (10) 

by taking the functional derivative with respect to the two-point function G, 

^ab{x, y) = -. (11) 

dGab{x,y) 

In practice, deriving S is not possible without making approximations since r2 consists of an 
infinite number of 2PI diagrams constructed from the bare vertex as well as the full 2-point 
function G. In Eq. (flOl) . Gq^ is the inverse classical propagator obtained from the classical 
action 5'[0], which defines the model to be considered, as 

igoa.(^,l/) = ,, ffit^, V (12) 
6(l)a{x)6(f)b{y} 

In the following we choose the 1/J\f expansion to next-to-leading order (NLO) where r2 is 
given as 

r2[0,G] = rLO[G] + rfO[0,G]. (i3) 

The leading order (LO) contribution to r2 consists of a 0-independent two-loop graph [fT3l 

r^°[G] = -^Jjr{Gix,x)}tr{Gix,x)} (14) 
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r^Lo [G] = 



r/Lo [G] = 








+ ... 



Figure 1. (Color online) The diagrams contributing, in LO and NLO of the 1/A/^-expansion, 
to the 2PI part r2 [0, G] of the 2PI effective action. At each vertex, it is summed over double 
field indices and integrated over double space-time variables. The wiggly lines represent the 
bare potential, the blue solid lines the full propagator G. 



which is shown in Figure da). The NLO contributions, Fig. Kb), can be summed up 
analytically, yielding 



with 



r^LOfG] = -Tr{lnB[G']}, 



B{x, y- G) = 5{x -y) + j^tr{Gix, y)G{y, x)}. 



Inserting Eqs. (fT3l)-(fT6l) into (fTTI) and using Eqs. (|7]) and ([8]) one obtains 



g[F{x,x) + -tT{F{x,x)}], 



g\I {x,y)F{x,y) - -F{x,y)p{x,y) 



^''{x,y) 

'LP{x,y) = - g^F{x,y)F{x,y) + I^{x,y)p{x,y) 
The functions I^'^ are obtained in analogy to Eq. ([8]) from the function 

H^^y) = ji= i^{G{x,y)G{y,x)} -i jj{x,z;G)tT{G{z,y)G{y,z)} 
which results from taking the derivative of r^'"^ with respect to G and are given as 



(15) 
(16) 

(17) 
(18) 
(19) 

(20) 



F{,x,yf - ^p{x,yf 
yo 



xo / 1 

2 ^ ^( ^ 



zS) 



F{x,z)\F{z,yY--p{z,yy 



2 1 I {x,z)Fab{z,y)pab{z,y) 



P{x,y) 



^^ 



Fab{x,y)pab{x,y) 



Xo 



F{x,z)Fab{z,y)pab{z,y) 



z,yo 



(21) 
(22) 



where = Fab Fab etc., and where a decomposition of / similar to that in Eq. ^ has been 
applied. 

This approximation has been shown to describe thermalization of a one-dimensional 
Bose gas starting far from thermal equilibrium lfT3l[T4l| as well as for relativistic field theories 
im [38l m Sni [lOl im. it includes off-shell scattering processes which allow the system to 
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exhibit damping, even in one spatial dimension. We point out, that the dynamic equations sum 
over the full past time evolution of the correlation functions and therefore include memory. 
Hence, no derivative expansion has been applied yet, which will constitute the main step to 
arrive at transport equations. 



3. Transport equations 



In this section we derive the transport equations in leading order (LO) and next-to-leading 
order (NLO) of a gradient expansion. We rewrite the exact dynamic equations given in the 
last section in order to get evolution equations with respect to centre and relative coordinates, 
providing a starting point for a gradient expansion and transformation to Wigner space. For 
simplicity we consider a system with spatially homogeneous initial conditions for the two- 
point function: 

F,fe(0, X, 0, y) = Fo^abi^, y) = 5,Jo(x - y). (23) 
The initial values of the spectral function are fixed by the commutation relations, cf. Eq. ([4]): 

Po,afe(x, y) = -ia2,abS{x - y). (24) 

To meet the requirements of the Fourier transformed equations in Wigner space on the 
one hand and to benefit from the spatial homogeneity of the system on the other hand, we 
rewrite the equations in terms of relative and centre coordinates 

X + y 



X 



x-y <^ X = X + s/2, y = X - s/2 



(25) 



where the two-point functions are to be transformed as, e.g., F{x, y) ^ F(X, s), z) i— > 
S^(X + s'/2, s - s'), F{z, y) ^-> F{X + {s' - s) /2, s'). We additionally introduce 

s' = z-y (26) 

which later on serves as the integration variable in the memory integrals. 

Using these definitions in Eqs. ([5]), Q, and adding these equations to the corresponding 
equations for the time derivatives with respect to y^, one obtains differential equations with 
respect to the centre time coordinate X^. 



dxoF{X, s) = io-2 



+ 



Hib{X + -) 



Hib{X - -) 



sW(X + ^)-SW(X-|)]f(X,s) 



e{Xo + 4 - 
- C'ix + 



2 ' 



S^(X + -,.-/)F(X 



s — s 



,s-s')j:^x 



-f{x+'-,s-s')j:^{x- 



2 

s — s' 



(27) 
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+ eW(x + ^)-s(°)(x-|mx,.) 



2 



p(X + -,s-s')S^(X- 



2' ' ' 2 



2 



- G^iX +'-,s- s')T.^{X - s')]]. (28) 

Note that, introducing retarded and advanced Greens functions and self-energies, 

G^{x, y) = e{xo - yo)pix, y), y) = -e{yo - xo)p(x, y), (29) 

S^(a;, y) = ^(xo - yo)^'ix, y), y) = -e{yo - xo)S''(x, y), (30) 

allowed us to send the integration limits of Sq in Eqs. (l27l) . (l28l) to ±oo. While the above 
equations have been obtained by adding the equations for dxoF{x, y) and dygF{x, y), etc., a 
second set of equations for the derivatives of s) and p{X, s) with respect to sq results 
when subtracting the respective expressions. These equations relate the two-time Greens 



functions to the single-time ones and are provided and discussed further in Appendix A 



3.1. Approximations 

With the aim to derive transport equations, we apply the following approximations: 

(/) The ^-function is neglected in the evolution equations for F, Eqs. (ITTI) . (lA.ll) . taking 
into account that the correlations disappear for large relative times. This corresponds to send- 
ing the initial time to to the infinite past. Note that, since an interacting system could have 
reached equilibrium at any finite time, transport equations are initialized by specifying F and 
p at a finite time using the equations with to ^ — oo as approximate description ll32l . 

(//) We apply a gradient expansion with respect to the centre coordinates X. 



3.2. Gradient expansion 



In leading order of the gradient expansion, taking into account the homogeneous initial 
conditions, one obtains from Eqs. (ITTI) and (l28l) the evolution equations 



dx,F{X,s) 
dxop{X, s) 



1(72 

0. 



F{X, s - s')J:p{X, s) - p{X, s - s')S^(X, s') 



, (31) 
(32) 
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Up to next-to-leading order in the gradient expansion these equations receive additional 
corrections as follows: 

dx,F{X, s) = LO + ia2^[sodx,T.^''\X)]F{X, s) 

^ ^ ^[49xoS+(X, s - s')]F{X, s') - S+(X, s'){s'^dx,F{X, s - s')] 



- Wx^iX, s - s')]E^{X, s') + s')[s',dx,^^{X, s - s')] 



(33) 



9xoP(X, s) = ia2|K9xoS(o)(X)]p(X, s) 

[49xoS+(X, s - s'MX, s') - s')[s',dxoP{X, s - s')] 

- [s',dxy{X, s - s')]J:f{X, s') + p+{X, s')[s',dx,^''{X, s - s')] 



(34) 



where 



p+ iX,s) = {X,s) + (X, s) , (35) 
S+(X,s) = S^(X,s) + S^(X,s) (36) 

and LO denotes the leading-order terms given by the right-hand side of Eq. (|3T1) . 
3.3. Transformation to Wigner space 

The transformation to Wigner space involves a Fourier transformation with respect to the 
relative coordinate s 



F{X,p) = I e'P'F{X,s), t^{X,p) = I e'P"S^(X,s), 

J s J s 



(37) 



etc., where ps = pqSq — p ■ s. 

For a spatially homogeneous system one finds that the diagonal matrix elements of F are 
purely real, while the off-diagonal matrix elements are purely imaginary. Furthermore, we 
find F'^{X,p) = F*(X,p) = F(X, —p). For p, the diagonal matrix elements are imaginary 
and the off-diagonal elements real. As compared with F, the function p changes sign under 
the transposition, p(X,p)'^ = — p*(X,p) = —p{X,—p). 

We apply the transformation to the previously derived equations of motion and have to 
take into account the integration limits when interchanging the derivative with respect to the 
centre time with the relative time integration. We thus introduce the derivative operator 

p2Xq n 

Vx,[F{X,p)]= / dso &'s(^^'dx,F{X,s) 



J-2Xo J 

dx,nx,p)-2 



F{X, so = 2Xo,p) 



+ e 



-2ipoXo 



F{X, so = -2Xo,p) 



(38) 
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Assuming lim^g^oo i^(so) = = lim^o^oo p(so) yields limxo-.oo ^^Xo = ^Xo- The leading- 
order transport equations result as 

Vx,[F{X,p)] = ia2{FiX,p)t'>{X,p) - p{X,p)t^ {X,p)] (39) 

Vxo[piX,p)] = 0. (40) 

The contributions from the integration limits account for the fact that the correlation functions 
are initialized at some initial time xq = yo = 0. These contributions can be removed for 
sufficiently late times as the time correlations are expected to vanish for sufficiently large 
relative times. 

To next-to-leading order in the gradient expansion, the transport equations read 



VxolF] = LO + a,\ [dxo^^'\X)]d,MX,p) 
1 



{ S+ (X, p) , F (X, p) }o + { (X, p) , p+ (X, p) }c 



2 

1^Xo[~P] = aJ[dx,^^'\X)]d,,piX,p) 



(41) 



- 1 (X, p) , p(X, p) }o + {E^(X, p) , p+ (X, p) }o' 
where we introduced the Poisson brackets with respect to Xq and po: 



(42) 



dA dB dA dB 
4. Comparison of dynamical and transport equations 

In the following we compare the transport equations derived in the preceding section, in 
leading order, Eq. (|39| ). and next-to-leading order, Eq. (HTI) . of the gradient expansion, with 
the equations of motion ([5]). This is achieved in two steps. First, we calculate F and p using 
equations ([5]), © in the NLO 1 /AT-approximation. We then take the results of this calculation 
and, after a Fourier transformation to Wigner space, we calculate both the left-hand and the 
right-hand sides of equations (|39l ) and (|4TI) . The left-hand side yields the time derivatives of 
the solutions of the 2P1 dynamic equations with respect to the centre time coordinate, while 
the right-hand side generates the corresponding approximative derivatives in leading order 
and next-to-leading order of the gradient expansion. 



4.1. Physical setup 

We consider a dilute homogeneous one-dimensional (D = 1) gas of spinless sodium atoms 
with mass m = 22.99 u (^'^Na) confined in a periodic box of length L = Nsttg ~ 43yU.m. 
The spacing of the numerical grid is = 1.33 pm, and the number of lattice points is 
Ns = 32. This results in discrete momentum mode^il p = {2Ns/ L) sm{np7r / L), with 

I Since the symbols x, p etc. so far have been defined as x = (xo,x), p = {po, p), we retain the boldface 
notation for space-like quantities although we consider the one-dimensional case. 
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rip = —Na/2...Ns/2. Using different initial line densities rii between 10^ and 10^ 
the total particle number varies in the range from ~ 4 to ~ 400 particles. 

For our investigations we consider different interaction strengths characterized by the 
parameter 7 which is related to the ID coupling constant as 

9 = (44) 

m 

We vary 7 in three steps from weak to strong interactions, 7 = 1.5 ■ 10~^, 0.15, 15, thereby 
keeping the quantity 'yn\ = mgrii fixed by choosing the line densities given above. In this 
way, the effect of quantum fluctuations is increased with 7 while the dynamics in the classical 
statistical approximation remains unchanged [14j. 

The spatial homogeneity allows for a Fourier transformation of the dynamic equations 
(in (6]) with respect to the spatial relative coordinate s to relative-momentum space: 

9^oF(xo,yo,p) = io-2|M(xo,p)F(xo,?/o,p) 
xo yo 
+ j dzoEP{xo,zo,p)F{zo,yo,p) - J dzo S^(xo, ^0, p)p(^o, l/o, p)|,(45) 


dxop{xo,yo,p) = i(T2|M(xo,p)p(xo,?/o,p) 

+ j dzoJ:P{xo,zo,p)p{zo,yo,p)^, (46) 

yo 

where the generalized energy matrix is given by 

2 r 

M{xo, P) = + I [tr{F(xo, xo, k)}l + 2F(xo, xo, k)], (47) 

according to (flTl) and [|T3ll . Here, Jj^ = (27r)"^ J dk denotes the one-dimensional momentum- 
integral, and the self-energy contributions are obtained, in NLO of the 1/J\f expansion, from 
the Fourier transforms of Eqs. (fTSi) to (|22l) which are provided in Appendix^ cf. Eqs. (|B.1I) - 
(Ell). 

With initial values for F(0, 0, p) and p(0, 0, p), the above coupled system of integro- 
differential equations yields the time evolution of the correlation functions. To fix F(0, 0, p), 
Eq. (l23l) . we chose a Gaussian momentum distribution peaked around p = 0, with a width of 

cr = 6.5 ■ 10^ m"^: /o(p) = "".(O, p) + ^, where 

n(0,p) = ^e-P'/"^'. (48) 
VTTcr 

The Bose commutation relations fix, as in Eq. (l24l) . the initial values for the spectral 
function: p(xo,a;o,p) = — io"2. 

4.2. Numerical solution of the dynamic equations 

The dynamic equations fl31) . (l46l) were solved using the techniques described in Refs. lfT3l[T4l| . 
From their solutions, the two-point functions in Wigner space were obtained by a discrete 
Fourier transform along the relative time direction sq. 
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Figure 2. The statistical and spectral two-point functions and self-energies for two different 
momentum modes as a function of the relative and centre-of-time coordinates sg Xq, 
respectively. Left column: momentum mode np = 0, right column; np = 8. From top 
to bottom: F, p, E^, and S''. The figures illustrate that all correlation functions vanish for 
sufficiently large relative times. See the main text for the parameters chosen. 



We expect the transport equations to reproduce, to a good approximation, the full 
dynamical evolution, if the two-point functions vary slowly with respect to the centre of time 
coordinate Xq as compared to the change with the relative time coordinate Sq. Figure [2] shows 
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Figure 3. (color online) Spectral function p (green solid lines), and statistical function F, 
normalized to the occupation number of the momentum mode p (blue dashed lines), as 
functions of po, for different times Xq and different momentum modes The difference 
between the respective spectral and normalized statistical functions decreases, indicating the 
emergence of a fluctuation dissipation relation. The good correspondence for the zero mode 
np = at the initial time is related to the choice of the initial condition. 



F(Xo, So, p) for two different momentum modes, together with the corresponding evolution 
of p, and S^. We see that for sufficiently late times the functions fall off to zero with 
increasing relative time Sq. The rapid decay of the self-energy as compared to F and p can 
be understood by noting that, in the next-to-leading order 1 /M expressions, see Eqs. (fTSi) and 
(fT9l) as well as (|2T)) and (l22l) . the correlators F and p enter to the third power. Compared 
to their oscillations along the relative time direction, we note only a weak dependence on 
the centre-of-time coordinate. The form of the correlators in the temporal plane shown in 
Fig.[2]is to be compared with that of the correlators for an ideal gas which show an undamped 
oscillatory dependence in the sq direction, with the frequency given by the free dispersion 
Pq = pV2m, see Eqs. (IA.14L (IA.16I) . 

As in the case of an ideal gas discussed in Appendix A. 2 the two-point functions F and 
p of an interacting system are expected to be related by a fluctuation dissipation relation, 

F{X,p) = ia2fn(p) + l)p{X,p), (49) 



2 

if the system is close to equilibrium [|4l|45]|. In Figure [3] we explicitly show that the function 
p2i {Xo, p) approaches Fn (Xq, p) / (n(Xo, p) + |) during the time evolution described by the 
2PI dynamic equations (|45] ), (|46l ), where n(Xo, p) = Fii{Xq, sq = 0, p) — |. 
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Figure 4. The momentum-mode occupation numbers n{t, p)/niL normalized by the total 
number of atoms in the box, niL = 853, as functions of time. The gas is in a far-from- 
equilibrium state initially, characterized by a Gaussian distribution ?i(0, p), Eq. ( |48] |, with 
width (7 = 1.3 • lO^m^^. It is weakly interacting, 7 = 1.5 • 10^^. Shown are the populations 
of the modes with p = 2Ns/ Lsnv{npTi /Ng), Up ^ (uppermost curve), 2,4, Ns/2 
(sequentially underneath), for Ng = 64, and one has n{t;—p) = n{t;p). A fast short- 
time dephasing period is followed by a long quasistationary drift to the final equilibrium 
distribution. Notice the double-logarithmic scale. 



4.3. 2PI dynamical versus transport equations 

We compare the time evolution as derived from the 2PI dynamic equations to their kinetic 
approximations following the procedure stated at the beginning of the section. The absolute 
values of the left-hand side (LHS), the right-hand side in leading order (LO) and the right-hand 
side in next-to-leading order (NLO) of Eqns. (l39ll4T)) are drawn in Figure [51 

Three generic time regimes are found O: Strong oscillations characteristic for 
early times, slow drifting at intermediate times, and a late-time approach to equilibrium 
characterized by vanishing time derivatives. These time regimes can already be seen in 
the corresponding time evolution of momentum-mode occupation numbers n(Xo, p) = 
Fii(Xo, So = 0, p) — I shown in Fig. |4] for a corresponding system with Ng = 64. Cf. Ref. 
lfT4l for a more detailed discussion of the numerical evaluation of the 2PI dynamic equations. 

The different regimes can be understood as follows: Close to the initial time, the 
strong oscillations of F{Xo,so,p) and p{Xo,so,p) in the centre time Xq are due to the 
finite integration limits Sq = ±2Xq in the sq direction. Consider the analytic solutions 
for the interaction free case, Eqs. (IA.14I) . (IA.16I) . Here, the oscillation frequency scales 
oc while the amplitude is constant, such that correlations do not decay. Consequently, 
the corresponding quantities in Wigner space, F{Xq,pq, p) and p(Xo,po) p)^ show the same 
oscillations if we restrict the Fourier transformation to the finite range —2Xo < sq < 2Xq. 
Besides these oscillations the amplitude of the correlation functions, in the interacting case, 
also decays with increasing sq. 
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Figure 5. (Color online) Time-derivative of the statistical two-point function as a function 
of centre time for two different momentum modes, rip = 4 (left panel) and = 8 
(right panel): Comparison of results from 2PI dynamic equations and those from their kinetic 
approximation. The (red) solid lines correspond to the left-hand side (LHS), the (blue) 
dashed(-dotted) lines to the the right-hand side in leading order (LO) (next-to-leading order, 
NLO) of the gradient expansion, Eqns. ( l39l l and dTTT i. respectively. In the left column of each 
panel, the frequency po is that of the peak of the spectral function, cf. Fig. [3]('on-sheir), in 
the respective right columns, po has been chosen two half- widths away from it ('off-shell'). 
From top to bottom, the line density rii is rescaled as well as the interaction parameter 7, so 
that mgni = jrif is kept fixed: Top row: case Ani ~ lO^m^^, 7 = 1.5 • 10^'^; middle row: 
case B ?ii = 10^m~^, 7 = 0.15; bottom row: case C ni = 10'''m^^, 7 = 15. 



The intermediate drifting regime is reached (for 7 = 1.5 ■ 10"^ at Xq ~ 0.003 s) when 
the contributions of the integration limits can be neglected. Non-vanishing derivatives of F 
and p with respect to Xq now solely result from the evolution of the correlations F{Xq, sq, p) 
and p{Xo, sq, p) in Xq. As can be seen in Figure [2l the change with Xq is slow compared to 
that with Sq. Finally, equilibrium is approached when all time derivatives vanish. 

We consider in more detail the cases corresponding to (A) weak, (B) moderate, and (C) 
strong effective interactions between the atoms: 

Case A corresponds to weak interactions, 7 = 1.5 ■ 10^^, i.e., to an essentially classical 
statistical evolution [[T4|. The evolution is expected to be well described by Boltzmann- 
type equations after the initial oscillations have damped out. The initial line density is set 
to rii = 10*^, which corresponds to ~ 420 particles. Figure[5K (top row) shows the left-hand 
side (LHS, solid line) as well as the right-hand side of the transport equations in leading order 
(LO, dashed line) and next-to-leading order (NLO, dashed-dotted line) of the gradient expan- 
sion for two different momentum modes Up = A (left panel) and Up = 8 (right panel). In 
the respective left columns, the frequency po is that of the peak of the spectral function, cf. 
Fig- [3] ('on-shell'), in the right columns, po has been chosen two half-widths away from it 
('off-shell'). We find that the 2P1 dynamic and the kinetic equations in general give the same 
results only as soon as the occupation numbers do not change any longer at all. 
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Figure 6. (Color online) The figure shows the time dependent decay constant F as defined 
in Eq. dSTT i. where the time derivative dF /dX^ is replaced by the left-hand side (LHS, solid 
Une), the right-hand side in leading order (LO, dotted line) and in next-to-leading order (NLO, 
dashed line) of the gradient expansion, Eqns. ( [391I4TI ). for the on-shell and the off-shell case 
as a function of the centre time. From top to bottom, the line density rii is rescaled as well as 
the interaction parameter 7 with = const., compare Fig. |5] The three rows correspond 
to cases A-C as described in the text. Colors indicate the momentum mode: ?ip = 4 (red). 
Up = % (blue). Up = 12 (green). Non-constant (and negative) values for F result from the fact 
that the time regime of exponential approach to equilibrium is not yet reached. 



The next-to-leading order contributions depend on the derivatives of the correlation 
functions. Compared to the interaction free case (IA.15I) . the correlation functions in Wigner 
representation become smooth functions of po as can be seen in Figure [3l This results in 
continuous po-derivatives, and therefore in non-vanishing contributions in next-to-leading 
order. However, as shown by the red dashed line, no significant contributions to the transport 
equations are found after the decay of the initial oscillations. 

Case B We increase the interaction parameter to 7 = 0.15 while the initial line density 
is decreased to rii = 10®, corresponding to a total number of ~ 40 particles. Note that in 
this way g oc 7/21 increases by a factor of 10 such that quantum statistical correlations grow 
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in importance, see Ref. lfT4ll . Our results are shown in the second row of Fig. [51 We find 
equilibrium to be reached faster, in particular for the higher momentum modes, while the 
intermediate drifting regime observed in case A is reduced. For the lower momentum modes 
we get similar results as in case A while for the higher momentum modes there is an essential 
difference to the preceding case: correspondence between 2PI and transport equations is not 
reached until equilibration occurs. 



Case C. We finally choose strong coupling, 7 = 15, and decrease the line density to 
rii = lO^m"^ corresponding to a total number of 4 particles. Our results are shown in the 
third line of Fig. [51 With this, we find qualitatively similar results as in the preceding cases 
and a faster approach to an equilibrium configuration. 

To study in more detail the late-time behaviour, we assume that, at late times, the 
statistical correlation function decays exponentially to its equilibrium value with a decay 
constant T{p), 

F{Xo,p) ^ F{Xo = oo,p) + AF(p)e-r(P)^°, (50) 

where AF{p) is some constant independent of Xq. In order to be independent of the 
equilibrium value, we plot, in Fig.[6l the time-dependent expression 

dF[XQ,p)/dXQ 

where dF / dX^ is given by the LHS, as well as by the LO and NLO expressions on the RHS of 
Eq. ([391 ). for three different momentum modes. The graphs correspond to the same choice of 
parameters and momenta as in Fig. [51 We focus the range of times Xq to those where T{Xq, p) 
is settling to a constant, indicating the emergence of an evolution according to kinetic theory. 
The top row of graphs in Fig. [6l shows that the different mode evolutions are settling to an 
exponential decay at times between 0.07 and 0.09 seconds, i.e., when compared with Fig.[4l 
only when there are almost no changes seen any more in the occupation number evolution. 
Hence, despite the fact, that a fluctuation dissipation relation is established almost an order 
of magnitude in time earlier, the kinetic approximation becomes strictly valid only at very 
late times. Fig. [6lalso shows that it is in general not sufficient to take into account the LO 
approximation in the gradient expansion only. 

Our findings become even more pronounced in cases B and C, when the system becomes 
more strongly correlated. For the largest couplings (case C) even to NLO in the gradient 
expansion the non-Markovian results of the 2PI dynamics are recovered. 

We finally study the dependence of the decay constants V on the density rii of particles 
along the one-dimensional system. Fig. [71 shows F, extracted at times tkin indicated in the 
inset, for five different densities Ui. Error bars indicate the variation over the different 
momentum modes. The times tkin have been chosen such that V(t) is found to remain 
approximately constant for times t > tkin- We find an approximately linear dependence 
of r on Til which indicates that the source of damping is rather an off-shell two-body than a 
three-body scattering effect. 
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Figure 7. (Color online) The decay constant F as defined in Eq. (ISTT i as a function of the 
particle density 7ii, extracted from the solution of the 2PI dynamic equations at times f^in as 
indicated with (blue) diamonds in the inset figure, after which the evolution can be described 
to a good approximation by an exponential decay. In the inset, the corresponding times idrift 
at which the slow drift of the occupation number sets in are indicated with (green) circles. The 
error bars indicate the variation of F, ikin, and tdrift over the different momentum modes. The 
red line is a linear fit. 



The inset of Fig. |7] also shows the times at which the fluctuation dissipation relation 
between the statistical and spectral correlation functions starts to hold. The difference between 
this time and the respective tkin increases with growing ni and therefore with increasing 
effective interaction strength. 

5. Conclusions and outlook 

We have studied the far-from-equilibrium dynamics of an ultracold, one-dimensional Bose gas 
and focused on the comparison between a fully dynamical approach on the basis of the 2PI 
effective action and the corresponding kinetic approximation in the form of Boltzmann-type 
transport equations. 

The 2PI effective action which allowed to derive dynamic equations by use of Hamilton's 
principle was considered in next-to-leading order of an expansion in inverse powers l/M 
of the number of field components. This approximation does not rely on a small coupling 
constant as expansion parameter and therefore is non-perturbative. Dynamic equations 
derived from the 2PI effective action automatically preserve essential conserved quantities 
as the total particle number and the total energy of the system. 

The transport equations were derived from the 2PI dynamic equations for the two-point 
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correlation functions by means of a gradient expansion with respect to centre coordinates 
and a subsequent Fourier transformation with respect to relative coordinates, i.e., a Wigner 
transformation. Furthermore, the details of the initial state are neglected by sending the initial 
time to minus infinity. This approximation relies on the fact that the time evolution of the 
correlation functions with respect to centre time occurs relatively slowly as compared to the 
oscillations with respect to relative time. 

We considered a homogeneous one-dimensional system of spinless bosonic atoms with 
an initially Gaussian momentum mode distribution. Both, weakly and strongly interacting 
systems were considered, such that the significance of quantum fluctuations to the time 
evolution could be studied. A comparison of the statistical and spectral correlation functions 
allowed to determine the time scale at which they are connected by a fluctuation dissipation 
relation. This timescale was shown to coincide with the time scale at which the transport 
equations set in to represent a valid description of the dynamical evolution. 

For weak couplings we observe good correspondence of dynamic and transport equations 
after an initial period of oscillations. However, off-shell effects are not covered by transport 
equations before equilibration occurs. Increasing the dimensionless interaction strength 
parameter 7 led to significant differences between the 2PI dynamic and the transport 
equations. Our results show that the kinetic description of the time evolution in terms of 
an exponential decay with decay rate V generically sets in to be valid only at very large times, 
when no essential change in the momentum profile of the system occurs any more. Moreover, 
they indicate that the late-time evolution is predominantly due to two-body off-shell scattering 
effects. 
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Appendix A. Transport equations in Wigner relative time directions 

Appendix A.I. Transport equations 

In Section |3] we have derived transport equations for the evolution along the Wigner centre 
time Xq by adding the dynamic equations governing y) and dy^F{x, y), respectively, 

cf. Eqs. (ITtI) and (l28l) . A second set of equations for the derivatives of F{X, s) and p{X, s) 
with respect to Sq results when subtracting the respective expressions: 




+ SW(X + |) + SW(X-^)]F(X,.) 
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+ S^(X + -,.-.')G'^(X- 2 



,^0 



(A.l) 



1(T2 



2' ' 2' 

+ sW(X + ^) + S(°)(X-|)]p(X,s) 

[S^(X + ^,.-s')p(X-^,.') 



p(X + -,.-s')S-^(X 



,^0 



s — s 



+ S^(X + -,.-s')G^(X 2 



G«(X + |,.-.')S''(X-^,.0] 



(A.2) 



In analogy to the derivation of the gradient expansion with respect to centre coordinates X in 
Section [J2l one obtains LO transport equations in the relative time sq, 



2m ^ ' 



F{X,s) 



\ s - s')F{X, s') + p+(X, s - s')S^(X, (A.3) 



dsop{X, s) = ids 



2m 



^ [S+(X, s - s')p(X, s') + p+(X, s - s')] |, (A.4) 



and their NLO corrections as 

d.,,F{X,s)=LO 



+ 



1(72 



[s'odx,mX, s - s')]F{X, s') - s')[s'odx,F{X, s - s') 



+ [s',dxoP{X, s - s')]J:^{X, s') - p{X, s')[s',dxo^^{X, s - s')]\ , (A.5) 

ds,p{X,s)=LO. (A.6) 

Here, LO denotes the respective right hand side of Eqs. (IA.3I ) and (IA.4I) . The functions p"*" 
and E+ are defined in Eqs. (l35l) and (l36l) . respectively. 

Proceeding with the transformation to Wigner space as in Sect. 13. 31 we obtain the LO, 



Vp,[F{X,p)]= ia2 



P 

2m 



sW(x)1f(x,p) 



(A.7) 



2PI nonequilibrium versus transport equations 



21 



as well as NLO equations: 



+ 2 



(^2 



t+{X,p)p{X,p) + p+{X,p)t'^{X,p) 



V,,[F] = LO -^[{E''{X,p),F{X,p)}o - {S^(X,p),p(X,p)}o 
Ppjp] = LO, 
where the operator Vp^ is defined in analogy to Eq. (|38] ): 



(A.8) 

(A.9) 
(A. 10) 



2Xo 



Vp,[F{X,p)]= / dso d'se^^'ds,F{X,s) 

J-2Xo J 



2Xo 

so=-2Xo 



(A. 11) 



Appendix A.2. Noninter acting gas 

In order to obtain an estimate for the form of the correlation functions in the Xo-So-plane, 
we derive analytical solutions for the interaction free case. This also recovers the fluctuation 
dissipation relation between F and p. 

Since the previously made approximations only affect the interaction terms, the leading- 
order transport equations are equivalent to the 2PI dynamic equations as can be seen by 
setting E(X, = 0. For a homogeneous system, spatial derivatives with respect to the 
centre coordinates disappear. Moreover, the momentum distribution /(p) = J d^s e~'P'^/(s) 
is real-valued. From equations (IA.7I) . (IA.8I) . we obtain 

^2 r n 2Xo 



-ip^F{X,p) = ia2^F{X,p) - 
2m 

These relations, together with the initial conditions (|24l) imply that 

,2 



e'^»^°F(X,so,p) 
e'P«^V(X,so,p) 



so=-2Xo 
2Xo 

so=-2Xo 



p{X, so,p) 
and, in frequency space, 

p(^,Po,p) 



, P 

-10-2 exp I l^CTsSo 



(A.12) 
(A.13) 

(A. 14) 



2m 



Po 



2m 



Po 



X 



{e{po)-9{-po))l + a2 



(A.15) 



where 1 denotes the 2x2 unit matrix. Since the equations (IA.12I) and (IA.13I) for F and p 
are identical, one obtains, with the initial momentum distribution /(p) = n(p) + 1/2, see 
Eq. (l23l) . the fluctuation dissipation relation (|49|) between F and p, such that 

2 

F(X,so,p) = /(p)exp {^^(^2So)- (A.16) 
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In this appendix, we provide the momentum-space self-energies T,^lf{xo,yo;p) in l-i-l 
dimensions which enter the dynamical equations (l45l) . (l46l) . From Eqs. (fTSi) to (|22l) . one 
obtains, for a homogeneous system, with (pa = 0, by Fourier transformation: 



where 



^abixo,yo,p) = ~g / I^ixo,yo,p- k)Fabixo,yo,k) 



Kbi^o,yo,p) = - 9 



~ ■^I''{xo,yo,P - k) pab{xo,yo,k) 



I''{xo,yo,p- k)Fabixo,yo,k) 
+ I^{xo,yo,p- k)pabixo,yo,k) 



(B.l) 



(B.2) 



I^{xo,yo,p) 



f/{F.(.„.,„.p-.)f.(.o....*) 

- \pab{xo,yo,p - k)pab{xo,yo, k) 



dzQ F{xo, Zo,p- k) [Fabizo, yo, k - k')Fab{zo, yo, k') 



- -^Pab{,zo, yo, k - k')pab{zo, 2/0, k') 



yo 



dzoI^{xo,zo,p - k) Fab{xo,yo, k - k')pab{xo,yo, k') 



(B.3) 



I''{xo,yQ,p) = 9 j^\Fab{xo,yo,p- k)pab{xo,yo,k) 

-If dzoPixo,Zo,p- k)Fab{zo,yo,k- k')pab{zo,yo,k')\, (B.4) 

Jk' Jyo ^ 

Here, = (27r)^^ J dk denotes the one-dimensional momentum integral. 
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